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1.0  SIMMARY 

The  objective  of  the  present  research  investigation  is  to  theo¬ 
retically  analyze  the  electromechanical  coupling  of  piezoelectric  zinc  sulfide 
films  on  silicon  substrates. 

During  the  past  twelve  months,  considerable  efforts  have  been  made  in 
analyzing  the  experimental  results  which  were  obtained  in  FY  1979.  In  order 
to  accomplish  this,  a  computer  program  was  written  to  calculate  surface- 
acoustic-wave  electromechanical  coupling  constants  for  layered  media.  The 
calculation  included  four  configurations  of  the  excitation  of  surface  acoustic 
waves  in  different  orientations  of  silicon  using  cubic  and  hexagonal  ZnS 
films.  The  calculated  results  showed  that  the  best  acoustic  coupling  (k2)  for 
ZnS  films  is  10  x  10"^  which  is  an  order  of  magnitude  smaller  than  ZnO  films 
on  silicon.  Also  included  in  this  report  is  the  calculation  of  ZnO  films  on 
fused  quartz  and  various  orientation  of  silicon. 

2.0  TECHNICAL  BACKGROUND 

In  1973,  Inaba,  Kajimura  and  Mikoshiba  theoreticaly  analyzed  ZnS 
films  on  fused  quartz2  and  showed  that  the  SAW  electromechanical  coupling 
constant  (k2)  was  7%  which  was  three  times  larger  than  LiNbOj.  The  same 
authors  measured  the  coupling  efficiency  of  ZnS  films  on  boro- si  11 cate  glass 
at  frequencies  ranging  from  70  to  90  MHz.  The  maximum  film  thickness  tested 
was  0. lh/X  where  \  is  the  SAW  wavelength.  The  best  conversion  efficiency 
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for  a  50  finger- pair  transducer  was  20  dB  with  electrodes  at  the  interface  and 
no  ground  plane. 

In  FY  1979,  ZnS  films  were  magnetron  sputtered  on  7059  glass  as  well 
as  various  orientations  of  silicon  substrates  at  Rockwell.  It  was  found  that 
the  best  insertion  loss  obtained  was  20  dB  with  30  finger-pairs  of  interdigital 
transducers  on  the  surface  and  a  floating  ground  plane  at  the  interface.  This 
corresponds  to  a  k2  of  10  x  10"*.  Since  the  dispersion  data  of  ZnS  on  glass 
or  silicon  were  not  available  we  decided  to  put  the  remaining  effort  into 
analyzing  and  calculating  such  layered  media  dispersion  curves. 

3.0  THEORETICAL  RESULTS 


3.1  Theoretical  Analysis 

The  Japanese  authors  in  reference  1  used  the  data  published  by 
Berl incourt  et  al2  erroneously.  The  correct  published  value  for  e^,  the 
piezoelectric  constant,  is  shown  in  Table  1. 


Table  1 


Elastic,  Piezoelectric  and  Dielectric  Constants  of  ZnS 


Japanese  Value 

Rockwel 1  Cal cul atl ons 

Cu  =  10.46  x  1010  M/m2 

cn  =  10.46  x  1010  N/m2 

C12  =  6.53  " 

C12  =  6.53  " 

C44  =  4.61  H 

C44  =  4.6  - 

e14  *  0.417  (C/m2) 

e14  =  0.14  (C/m2) 

ell/co  =  8,32 

ell/£o  *  8,32 

p  *  4083  kg/m3 

p  *  4083  kg/m3 
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For  cubic  ZnS  crystals,  the  non-vanishing  matrix  for  elastic, 
piezoelectric  and  dielectric  constants  are  given  below 


El  astic 

C11 

C12 

C12 

C12 

C11 

C12 

Cij  =  C12 

C12 

C11 

•  •  •  C44  •  • 

.  .  .  .  C44 


Piezoelectric 


e14  * 

e14  • 

e14 


Dielectric 


G11  * 

ell  • 
G11 


Since  the  sputtered  ZnS  films  on  glass  are  generally  hexagonal  in 
symmetry,  a  rotation  of  the  matrix  elements  in  the  cubic  system  to  (111)- plane 
will  give  the  required  elastic,  piezoelectric  and  dielectric  constants.  The 
rotated  non-vanishing  matrix  elements  for  (111)  plane  are  given  in  Table  2* 
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Elastic 


Table  2 


Rotated  Elastic, 

Piezoelectric  and  Dielectric  Constants 

for  (111) 

Plane  of  ZnS 

C' 

hi 

C12 

C' 

h3  • 

• 

C12 

hi 

C13  * 

• 

Cij  =  C13 

C13 

C33  * 

•  • 

* 

• 

'  C44 

•  * 

°44 


"66 


Where 


U11 

C12 

q3 

C33 

C44 

C66 


13.11  x  10 
5.65  " 
4.77  " 
13.99  " 
2.847 
3.729 


10 


N/m2 


Piezoelectric 


•ij 


15 


'15 


e31  e31  e33 


Where 


-15 

?31 

?33 


-0.08  C/m2 
-0.081  " 
0.162  " 


Dielectric 


:11 


1j 


:11 


■m 


Where  /e^  =  8.32 

11  0 
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Consider  a  surface  wave  propagating  along  the  x-axis  of  the  layered 
media  shown  in  Fig.  1  with  X3  axis  normal  to  the  surface.  The  stress  equation 
is  given  as: 


(1) 


where  is  the  stress,  Uj  the  mechanical  displacement  and  p  the  density. 
The  electric  displacement  for  an  insulating  material  such  as  ZnS, 


(2) 


For  piezoelectric  material,  the  equations  of  motion  are  given  as: 


Tij  “  Cijkx.  \  1  "  ekij  ^ 
Di  =  eik£  Skt  +  eik  ^ 


Where  Sk1  =  1/2 ( 3Uk/3xA  +  aU^/Sx^)  are  components  of  strain,  =  3<i>/3xk  the 

electric  field,  the  elastic  constants  at  constant  electric  field,  e^j 

c 

the  piezoelectric  constants  and  the  dielectric  permittivity  tensor  at 
constant  strain. 


Substituting  Eqs.  (1)  and  (2)  into  (3)  and  (4)  to  give 


iiu.i  -  e  »  ,8VSV.  A 

aTp  ‘  7  ci jkt  Wj  1 +  3xj^  +  eki j  “  p  ~^T 


(5) 
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J. 


,  3U.  3U0  c  .2. 

_ ( !  4  21  \  _  p  ^  ”  T  _ 

ikt  axi  3xk'  Ik  ax^ax^ 


=  0 


i,  j,  k,  t  =  1,  2,  3 


(6) 


Assuning  a  particle  displacement  of  the  form 
ikbx,  ik(x,  -  vt) 

U .  =  a .  e  e  (7 ) 

J  J 

where  aj  are  unit  vectors  along  Xj  respectively;  wave  member  k  =  2"A,  b  the 
decay  constant  with  depth  and  v  the  phase  velocity;  and  the  electric  potential 

as 

ikbx,  i k( x,  -  vt) 

<j>  =  e  e  ( 3 ) 


Substituting  Eqs.  (7)  and  (8)  into  (S)  and  (6)  gives: 


"rll  -  pv2 

r12 

r13 

r14 

r-  — 

ai 

r12 

r22  -  pv2 

r23 

r24 

• 

a2 

r13 

r23 

r33  -  pv2 

r34 

a3 

r14 

L. 

r24 

r34 

r44_ 

_a4 

=  0 


(9) 
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where 


rll  =  c55  1)2  +  20 15  b  +  C11 

r22  =  c44  1)2  +  ^46  b  +  C66 

r33  =  C33  1)2  +  ^35  b  +  C55 

r12  =  C45  1)2  +  ^C14  +  c56^  b  +  c16 

r13  =  C35  b2  +  *c13  +  C55^  b  +  c15 

r23  =  C34  b2  +  ^c36  +  C45^  b  +  c56 

r14  =  ®35  f)2  +  {e15  +  e31)  b  +  ell 

r24  =  e34  b2  +  ^e14  +  e36^  b  +  e16 

r34  =  ®33  1)2  +  ^e13  +  ®35^  b  +  e15 

with  the  electric  term 

r44  =  -<e33  b2  +  2  £13  b  +  e^) 

The  secular  Eq.  (9)  is  applied  to  the  layer  and  substrate  regions  and 
the  matrix  of  coefficients  is  reduced  to  a  polynomial  in  b.  In  the  completely 
general  case  the  polynomial  is  of  3th  order  and  8  complex  roots,  bj  are  solved 
for.  For  the  substrate,  only  those  roots  which  cause  decay  to  zero  with 
decreasing  x3  are  retained.  This  results  in  12  solutions  to  the  equation  of 
motion  (3). 

The  total  solution  to  the  layered  problem  is  taken  as  a  sum  of  the  12 
partial  solutions. 


-  In)  .,kbnx3  .,k(xl  '  *« 
Cn  aJ  *  * 


(10) 
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♦  = 


z 

L»S 


(n) 


ikb  x- 
„  n  j 


ik(Xj  -  vt) 


n  =  1,2,3 


where  \  indicates  summation  over  root  nunbers,  n,  for  the  layer  (L)  and 
L.S 

substrate  (S). 


The  general  boundary  conditions  are: 


a)  Mechanical  Transverse 

1.  Continuity  of  transverse  displacement  at  the  interface 

Us  =  U  2  at  Xj  =  0 

?.  Continuity  of  transverse  shear  stress  at  interface 

T32  =  ^32  at  Kj  =  0 

3.  Vanishing  of  transverse  shear  stress  at  the  free  surface 

T^  =  0  at  Xy  h 

b.  Electrical 

4.  Continuity  of  the  normal  component  of  electrical  displacement  at 
i nterf ace 

°3  -  0  *  at  X3  »  0 
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5.  Continuity  of  potential  at  interface 


ks 


at  ru  =  0 


6.  Continuity  of  normal  component  of  electrical  displacement  at 
free  surface 


°3  =  k  eo  ♦*  at  x3  -  h 


(for  °3  =  '  eo  %  when  V  -) 


c.  Mechanical  Sagittal 

7.  Continuity  of  longitudinal  particle  displacement  at  interface 


U*  =  uj  at  Xj  =  0 


8.  Continuity  of  vertical  particle  displacement  at  interface 


U3S  -  55 


at  x,  =  0 


9.  Continuity  of  sagittal  shea^  stress  at  interface 


T31  =  ^1  at  X3  =  0 


10.  Continuity  of  vertical  compressional  stress  at  Interface 


^3  *  ^33 


at 


*3  * 
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11.  Vanishing  of  sagittal  shear  stress  at  free  surface 

T^l  =  0  at  Xj  =  h 

12.  Vanishing  of  vertical  compressional  stress  at  free  surface 

=  °  at  x3  =  h 

If  the  surface  or  interface  is  conducting  then  the  conditions  on  displacements 
D3  are  changed  to  reflect  a  zero  potential  at  the  appropriate  value  of  X3. 

Thus  the  above  12  boundary  conditions  provide  12  equations  in  the  12  unknowns, 
Cp.  The  matrix  elements  of  the  general  boundary  condition  determinant  are 
shown  in  Table  3. 

In  order  to  obtain  non-trivial  solutions  to  this  set  of  homogeneous 
equations,  a  12  x  12  determinant,  successive  values  of  v  are  chosen  until  the 
boundary  condition  determinant  is  equal  to  zero. 

A  computer  program  based  on  the  above  analysis  was  developed  to 
calculate  the  dispersion  characteristics  and  electromechanical  coupling 
constants  of  ZnO  and  ZnS  transducers  using  the  piezoelectric  data  of 
Berl Incourt  et  al . 

3.2  Zinc  Oxide  Films  on  Fused  Quartz  and  Silicon 

As  a  test  of  the  computer  program  and  for  comparison  of  performance 
with  ZnS  films,  ZnO  on  fused  quartz  and  (111)  silicon  substrates  was  analyzed. 
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The  results  for  ZnO  on  fused  quartz  are  shown  in  Fig.  2.  Figure  2a  shows  the 
velocity  dispersion  as  a  function  of  layer  thickness  normalized  to  wavelength 
over  the  range  0.001  to  5.  The  normalized  velocity  change  (AV/V)  correspond¬ 
ing  to  a  transducer  (electrodes  on  top,  ground  plane  at  interface)  shown  in 
Fig.  2b.  These  calculations  were  checked  against  previously  published  results 
of  Kino  and  Wagers  [3]  for  comparison.  The  coupling  peak  at  5%  is  approxi¬ 
mately  5  x  10~^. 

A  similar  calculation  was  then  made  for  ZnO  on  (111),  <110>  silicon 
and  the  results  are  shown  in  Fig.  3.  The  maximum  coupling  peak  again  was  at 
5%  however  somewhat  reduced  to  3  x  10"^. 

3.3  ZnS  Films  on  Fused  Quartz  and  Silicon 

After  determining  the  correctness  of  our  ZnO  calculations,  the 
analysis  of  ZnS  films  was  performed.  The  first  material  analyzed  was  ZnS  on 
fused  quartz.  Surface  wave  velocity  and  coupling  were  computed  for  all  four 
possible  transducer  configurations  as  shown  in  Fig.  4.  The  peak  in  coupling 
at  5%  was  much  less  than  the  previous  ZnO  result,  i.e.,  1  x  10"*.  This  is  not 
equal  to  LiNbO-i  and  more  than  an  order  of  magnitude  less  than  ZnO  on  fused 
quartz. 

ZnS  films  on  (001)  propagation  along  <100>  and  (111),  <ll0>  silicon 
were  also  analyzed.  The  results  are  shown  in  Fig.  5  and  6.  On  (001)  silicon 
the  coupling  at  5$  was  1  x  10"*  and  on  (111)  the  coupling  was  0. 5  x  10"*. 
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h/A 


h/A 


Fig.  2  ZnO  film  on  fused  quartz;  a)  Phase  velocity  and  b)  electro¬ 
mechanical  coupling  constant  k2 (2  AV/V)  as  a  function  of  h/A 


VELOCITY  (M/SEC) 


VELOCITY  (M/SEC) 
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h  /X 

y  _ 

Fig.  6  ZnS  film  on  (lll)-plane  of  SI,  propagation  vector  along  <110> 
direction;  a)  phase  velocity  and  b)  electromechanical  coupling 
constant  (2  a V/V)  as  a  function  of  h/X. 
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4.0  CONCLUSIONS 

Theoretical  calculations  were  performed  on  ZnO  and  ZnS  films  on 
silicon  and  fused  quartz  substrates.  ZnO  typically  showed  a  peak  in  coupling 
of  5  x  10"3  at  5%  thickness  normalized  to  wavelength.  ZnS  films  were  found  to 
have  considerably  less,  typically  1  x  10-4  at  the  same  film  thickness.  This 
result  directly  conflicts  with  previously  published  calculations  [1]. 

However,  comparing  the  ZnO  and  ZnS  (hexagonal)  piezoelectric  constants  in 
Table  4,  the  result  is  expected.  This  is  substantiated  by  our  experimental 
measurements  on  ZnS  film  coupling. 


Table  4 

Piezoelectric  Constants  of  ZnO  and  Rotated  ZnS 

ZnO  ZnS 

e-jj  =  1.321  =  0.162 

eji  =  0.573  e^j  =  -0.081 

ejg  =  -0.480  ®15  =  “0*0818 
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